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Abstract 

A new model for the exchange of two gluons between the virtual photon and 
the proton, in non-diffractive deeply inelastic electron-proton scattering, is 
developed and studied. This model relies on a perturbative calculation, pre- 
viously applied to diffraction, and a general result from Regge theory. As a 
first application of the model, we study corrections to the momentum trans- 
fer to the quark-anti-quark pair, at the photon-vertex. We find a significant 
enhancement of the cross-section at ~ Q 2 momentum transfers, and large 
negative corrections for small momentum transfers. The implication of this 
result for jet-distributions measured at HERA, is discussed. 

PACS: 13.60.-r 
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1 Introduction 



In an additive parton model description, the rapid growth of the electron-proton 
cross-section with the total energy s, is interpreted as a growth in the number 
of partons (mainly gluons), inside the proton. When the density of the partons 
is sufficiently large, it is generally expected that they will start to self interact 
(screening corrections), and that the quark probe, which measures the number of 
gluons, will interact with several gluons simultaneously (shadowing). Such correc- 
tions are expected to slow down the energy growth of the cross-section, in a way 
that the Froissart bound, a oc log 2 (s) |L|, is not violated. 

It has, however, been shown, for large enough virtualities of the photon, Q 2 > 
1 GeV 2 , where perturbation theory is applicable, that the structure function mea- 
sured at HERA, is described both by models that do not include multiple interac- 
tion effects 0, and by models where multiple interactions have a significant con- 
tribution (See e.g. B). The main theoretical uncertainty, is the non-perturbative 
nature of the proton. Different perturbative behaviour can describe the same struc- 
ture function, by making different non-perturbative assumptions for the proton. 

One way of settling this ambiguity, is to study more exclusive observables, e.g. 
properties in the final-state distribution, in addition to the structure function. A 
striking example of this is the observation of rapidity gap events. The appearance of 
large regions in rapidity with no hadrons in the final state, are generally considered 
to be a signal of the exchange of a colour singlet state between the proton and the 
probe. In a perturbative description, this is only possible by the exchange of, 
at least, two gluons. Thus, the rapidity gap events, indicate the significance of 
multiple interactions. 

It is nowadays, well established that, in many cases, multiple interactions have 
a significant effect also on the non-diffractive final states. This became first evi- 
dent, in proton-anti-proton scattering. It was observed that multiple interactions 
considerably increase the fluctuations in the final-state multiplicities of charged 
hadrons. Assuming universality for the hadronisation process, it was shown |3J] 
that multiple interactions are necessary for the description of high multiplicity dis- 
tributions, measured by the UA5 collaboration for proton-anti-proton scattering at 
540 GeV [|]. Direct measurements of the cross-section for double independent par- 
ton scatterings, have been made, from the appearance of two uncorrelated pairs of 
balanced high-p^ jets ||[7] (In j7|, the measurement was made on 3 jets + photon 
final-states.). 

The same phenomenon is observed in the case of electron-proton scattering at very 
low virtualities of the photon. The hadronic component of the quasi real photon, 
enables a description of multiple interactions, similar to the case of hadron-hadron 
scatterings. Again, one picks out two (or more) partons from each projectile, and 
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describe the multiple interaction process as simultaneous 2^2 scatterings ||. In 
this way, the inclusive distribution of each scattering becomes the same as for single 
scattering, or in other words, inclusive final-state distributions are not changed. In 
photon-proton scattering, a clear signal for multiple interactions were found in the 
form of the so called "Pedestal effect" with an enhanced low-p_i_ activity in the 
surrounding of high-p^ jets |J. 

In this article, the independent multiple scatterings, described above, is comple- 
mented with a model which relies on a perturbative treatment of the coupling of 
the photon, to two gluons inside the proton. This model results in a different kind 
of final-state, compared to the case of two independent scatterings. It does not nec- 
essarily lead to an increase in the jet activity in the photon direction, but instead, 
it leads to a change of the inclusive jet distributions. It is perhaps also more rel- 
evant for virtual photons, which should not be expected to have non-perturbative 
hadronic properties. 

In Section ^|, the basic ideas behind the model are described. One of the main 
ingredients is a perturbative treatment of the photon vertex which is relevant for, 
and has previously been applied to, diffractive events in electron-proton scattering 
TOfl . A general result from Regge theory (The AGK cutting rules [jlTJ] ) , is then used 



to relate this perturbative calculation, to the corresponding double gluon exchange 
corrections to non-diffractive DIS events. By following this strategy, the relevant 
formulas (Eqs. ([12]) and (p~3|) ) are derived in Section [| In Section f|, the four-gluon 
impact factor of the proton is determined, by phenomenological arguments. 

In Section [5], the model is used to make numerical estimates for the effects on the 
momentum transfer (pj_) between the proton and the quark- ant i- quark pair at the 
photon vertex. We find that the differential cross-section da/dp 2 ^, for p\ ~ Q 2 is 
increased with about 10%, for Q 2 = 10 GeV 2 and xb = 0.0001, relevant for the 
HERA region, and that this correction grows rapidly towards lower values of xb 
and Q 2 . It is however pointed out that the uncertainties in this estimate are large, 
in this first treatment. In the low-pj_ region, we find large negative corrections. 
This is interpreted as a signal for saturation in the p^-distribution. 

In the Appendix, we make a consistency check for the validity of the AGK-rules, 
in the sense they have been used for the model. This is however not a full QCD 
calculation. 

It should be mentioned, that the same kind of physics that is discussed here, can 
be treated inside a formalism that has been developed for radiation process in 
medium So far, this study has been concentrated on the effect of energy 



loss in medium, but recently JL3||, it was extended to also cover effects on transverse 



momenta distributions. Explicit studies in this respect for e-p scattering, has 
however not yet been done. 
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(a) (b) (c) 

Figure 1: (a) Normal DIS: A single gluon is exchanged between the quark-pair 
(above) and the proton (below). The proton breaks up into a high-mass hadronic 
state, (b) N on- diffract ive double gluon exchange: The hadronic state becomes 
more dense, (c) Diffractive event: A pomeron is exchanged. It couples to the 
quark-pair as a system of two gluons in colour singlet. 





Figure 3: Examples of two-ladder exchange, with different final states: (a) Both 
ladders are cut and the final-state distribution is doubly dense, compared to normal 
DIS. (b) One ladder is cut. This gives the same final-state distribution as the 
leading contribution. 

2 A model for perturbative double gluon exchange 

The details of the present model are outlined in the following two sections. Here, 
we discuss the physical problem in general, and describe the main assumption that 
are used in the model. 

In the region of low xb and moderate Q 2 , the leading log(s) (LLs) calculations 
of BFKL [[TJJ are expected to be valid, and gave in fact, a qualitatively correct 



prediction for the growth of the cross-section, towards lower x^-values, a oc x^ x . 
(It should be mentioned that also models based on DGLAP evolution give a good 
description of the strong rise of the structure function ||.) In this frame work, the 
leading contribution to the cross-section is given by the exchange of a colour singlet 
gluonic ladder between the proton and a quark loop at the photon vertex (See Fig. 
H). The rungs in the ladder appear as gluon emissions into the final-state. After 
the hadronisation process, one will then end up with a hadronic final-state with a 
very large mass, W 2 = Q 2 /xb- The BFKL prediction for the value of the exponent 
A = 4alog(2) « 0.55, for a = Na s /n = 0.2 with N = 3 colours, is observed to 
be too large compared to measurements of the structure function. Sub-leading 



corrections are expected to lower this value to around A = 0.3 — 0.4 (See e.g. |15 



and |I6|). For this cross-section, down to the next-to- leading order, we will use 
the notation a\ , since it corresponds to the process that the quark pair absorbs 
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a single gluon from the proton and has its contributions given from diagrams with 
two reggeised gluons coupled to the quark loop. 

The multiple interaction effects that we are interested in here, appear first at the 
sub-sub-leading level and are currently not fully calculated. For double gluon ex- 
change! (With the cross-section given from diagrams with four reggeised gluons 
coupled to the quark loop.) one can however expect contributions with two ex- 
changed ladders, which would lead to a faster energy growth than the leading 
contribution. We will use results from Regge theory, as a way to determine just 
these fast growing contributions at the sub-sub-leading level. 

In the language of Regge theory, the gluon ladder which appears in the leading 
BFKL contribution corresponds to the exchange of a pomeron, with the final-state 
cut going across it. In Regge theory, one can also discuss the possibility of the 
exchange of several pomerons. Two-pomeron exchange, with the final-state cut 
going in between the two pomerons, has been successfully used for the description 
of diffractive events in electron-proton scattering ||10|| . There, the pomeron and its 



coupling to the proton, was determined by phenomenological arguments, while, at 
the photon vertex, the pomeron was assumed to couple as two gluons in the colour 
singlet state. The leading part of the amplitude for pomeron exchange is expected 



to be purely imaginary. According to the Cutkosky rule |T7| , this means that it is, 
to leading order, given by its s- discontinuity, 

A diS « ilnuAdiff = ^disc s A diS 

The diffractive cross-section, which we will denote cq, is in this way determined 
from a triple discontinuity and the calculations are considerably simplified. 

In Regge theory, the two-pomeron exchange gives two other types of contributions, 
but with non- diffractive final-states. These are given by the possibilities that the 
cut goes through one or both of the pomerons. We will base the modelling of 
double gluon exchange and a higher order correction to single gluon exchange on 
these non-diffractive contributions of two-pomeron exchange. The contributions 
to double gluon exchange, the process that the quark pair absorbs two gluons, 
are given from diagrams like the one in Fig. ||]a with two cut pomerons. The 
contribution from diagrams with one cut pomeron (Fig. |3|b), correspond to the 
higher order correction to single gluon exchange. 

Throughout in this article, we denote the cross-section of single gluon exchange a\ 
and the cross-section of double gluon exchange with <ji- Single gluon exchange has 



2 Note that we are here discussing the underlying physical process of two exchanged gluons, 
however, the processes of double gluon exchange and the correction to single gluon exchange, 
have their contribution given from Feynman diagrams with four exchanged gluons. We use the 
notion 'single' and 'double' when discussing the physical process, and numbers e.g. 'two-', 'four-' 
for the discussion of Feynman diagrams. 
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two separate contributions, u\ = + a[ 4 \ where is the BFKL cross-section 
with two reggeised gluons coupled to the quark loop, while is the sub-sub- 
leading correction, with four reggeised gluons. (Some sub-leading corrections in 
the gluon ladder are also included by using the lower value of the exponent A.) 
The diffractive cross-section is denoted cro- For the sake of discussion, a tilde will 
be used for the cross-sections defined in Regge theory. Thus, we will denote the 
two-pomeron exchange cross-section with i cut pomerons with 6^. However, in 
the present approach, each of these will be identified with the respective untilded 
cross-section (With the same sub-script). 

For the non- diffractive contributions of two-pomeron exchange, i.e. a 2 and a^, it 
is not certain that the use of the Cutkosky rule gives the complete answer, since 
the amplitudes (The diagram on either side of the cut in Fig. |3| a and b.) can have 
real contributions, also at the level, in the LLs expansion, that we are interested 
in. We will resolve this uncertainty, by using a result in Regge theory. In an 
article by Abramovski, Gribov and Kancheli it was shown that the relative 



size of the contributions of two-pomeron exchange, but with different number of 
cut pomerons is given by 

a Q : a x : a 2 = 1 : -4 : 2. (1) 

The main ingredient in their proof was factorisation into one, approximately imag- 
inary, signature factor for each pomeron, and a single vertex function for the 
coupling of the two pomerons to the external particles. However, they used an 
underlying 3 -theory of scalar particles, in order to show that the vertex functions 
are not changed in the different cuttings. We will make the non-trivial assumption 
that the AGK rules are valid also in our case. This allows us to use the diffrac- 
tive calculation, where the Cutkosky rule is applicable, to determine the other two 
contributions and a 2 - 

As described above, the diffractive cross-section is given from diagrams with one 
pomeron on each side of the final-state cut. Each pomerons couples to the quark 
loop as two gluons in the colour singlet state and with the transverse momenta 
relatively well balanced. For the case of double gluon exchange a 2 (Fig. |3]a) 
one also has two gluons that couples to the quarks, on each side of the final- 
state cut, but they arise from two different pomerons. They are not necessarily 
in the colour singlet state, and the momenta can be rather unbalanced, due to 
perturbative QCD radiation. Finally, the higher order correction to single gluon 
exchange, cr^ (Fig. |3]b) is given from diagrams with one gluon on one side of the 
cut and three on the other. In this case, the transfered momentum and colour 
is given from that of the single gluon. In all the three cases, we will assume 
that each individual gluon, identified by its origin from the pomerons, has the 
same momentum distribution as determined in the diffractive case, from the triple 
discontinuity calculation. In this way, we can calculate the momentum transfer also 
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for the two non-diffractive contributions. The normalisation of each contribution 
is determined from the AGK cutting rules (and the diffractive calculation), so that 
the relative size of the inclusive cross-sections will satisfy Eq. ([!]). 

The applicability of the AGK rules in QCD was discussed for the case of quark- 
quark scattering in the report article of Gribov, Levin and Ryskin Jl9| and the 



case of photon-proton scattering was treated more recently, in a article by Bartels 
and Ryskin ||18|| . This latter treatment was based on contributions from triple 
discontinuities, as for the diffractive cross-section discussed above, but for arbitrary 
colour states. 

In the present model we are assuming that there are further important contribu- 
tions, apart from the triple discontinuities, which would arise from the real parts of 
the relevant amplitudes^. We will not here attempt to calculate these sub-leading 
contributions from first principles, but the usage of the AGK cutting rules suggest 
that these are necessary in our treatment. 

In the Appendix, we show, for the symmetric colour exchange part of a%, that the 
triple discontinuity result completely agrees with the AGK rules as they are used 
here. We are thus assuming that o<i has important contributions from the real part 
of the relevant amplitudes, but only in the case that the two exchanged gluons are 
in a anti-symmetric colour state. 

We end this section with a more intuitive discussion of the above results, in par- 
ticular, the interpretation of the higher order correction and its AGK factor 
relative to 02. Consider a model, e.g. the BFKL formalism, which approximates 
the cross-section with the single gluon exchange cross-section, a\ . The inclusive 
cross-section can, in this approximation, be described as the sum of all possibili- 
ties to exchange a single gluon between the quark pair and the proton. However, 
the possibility of the simultaneous exchange of a second gluon is not excluded, 
but it goes beyond the accuracy of the calculation. In the leading BFKL calcula- 
tion, the correction is formally small, since it has two extra factors of the strong 
coupling but no new factors of log(s). This approximation leads to an overesti- 
mate of the total cross-section because, roughly speaking, events with double gluon 
exchange contribute twice. On the other hand, the one gluon exchange approxima- 
tion underestimates the momentum transfer to the quark pair, simply because two 
exchanged gluons in general carry more momentum than one. When double gluon 
exchange corrections are included, as in the present model, we have two separate 
kinds of events: We have a new kind of events with two exchanged gluons, with 
the cross-section cr 2 , and we have the single gluon exchange events with a modified 
cross-section o\ = af^ 1 + a^f 1 . The role of is to correct for the double counting 
(in crj 2 ^) of the double gluon exchange events. Therefore, one can intuitively expect 



3 Note that we here in general use the word "amplitude" for the diagram on either side of the 
final-state cut. 
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the AGK relation a\ ' = —1o 2 to be valid. 

There are two important implications from this. First, that the total higher order 
correction is negative. We know this because, a 2 corresponds to the probability of a 
certain process, and must be positive, while the total correction, = + a 2 = 
—a 2 is negative. The other implication, that we can expect and will become more 
explicit in the treatment below, is that the differential cross-section at a certain 
finite momentum transfer is safe towards the exchange of a very soft gluon. If one 
of the two exchanged gluons is very soft, a 2 will get a contribution at a momentum 
transfer which is close to the momentum of the harder gluon. The correction <j[^ 
on the other will get half of its contribution at the softer momentum and half 
of it at the higher one. At the larger momentum transfer, the two contributions 
will, in the limit, cancel exactly, while there will be an extra contribution, which 
we usually do not have good control over, around the point of zero momentum 
transfer. In a perturbative treatment, it is therefore better to study corrections to 
the differential cross-section instead of the inclusive. 



3 The coupling of two gluons to the photon 

Let us start with a more detailed discussion of the amplitude which is relevant 
for the exchange of two gluons in the photon-proton scattering process. On the 
upper side, we have a virtual photon which breaks up into a quark-anti quark 
pair. On the lower side we have the incoming proton. Two gluons are exchanged 
between the proton and the quark pair, either as two separate exchanged gluons or 
they can combine into one gluon which is absorbed by the quark-pair. The proton 
transforms into a hadronic state with high mass, in the case of non-diffractive DIS, 
and a lower mass state in case of pomeron exchange. 

The two gluons can be in either of the colour states of a singlet, an anti-symmetric 
octet or a symmetric octet. Higher representations are not possible for the gluons 
in this process since the outgoing quark pair lies in the 3 <g> 3 = 1 © 8 colour space. 
We will now consider the different colour states separately and describe them in the 
frame work of the leading logs expansion (LLs). Denote the amplitude for singlet 
exchange Ai, for anti-symmetric octet exchange Ag a and for the symmetric octet 
Ass- The leading part of As a is known. It factorises into the amplitude for one- 
gluon exchange and, to 0(a^), a factor u(t)\og(s) where uj is the trajectory function 
of a reggeised gluon and t is the momentum transfer. In the partonic language, it 
corresponds to the exchange of a single gluon, and the factor u(t)log(s) corresponds 
to a virtual emission. In our case we will be interested in contributions that belong 
to the next order of the LLs expansion, that is, they lack any factors of log(s), but 
contain two factors of the strong coupling a s . 
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In the following, it will be useful to separate the amplitudes into their real and 
imaginary parts. The reason is that the imaginary parts are (for asymptotic s) 
easily calculable through the Cutkoski rule. We can then use the AGK cutting 
rules to also determine some contributions, due to the real parts of the amplitudes, 
which are important in the model. In Regge theory, it is expected that both of the 
amplitudes A\ and Ag s correspond to the exchange of even signature reggeons. The 
leading part of a reggeon with even signature is expected to be purely imaginary, 
thus we can assume Ai = ilnuAi and A$ s = ilmAss- The leading part of As a is 
real (as mentioned above) while the next-to leading contribution has both a real 
and an imaginary part. In our case, when the two gluon exchange amplitude is 
multiplied with its conjugate, the contributions A$ s Al s and Im^4 8a Im^lg a 

can be calculated in the asymptotic limit, using the Cutkosky relation, while the 
relevant contributions from Re^^^Re^g^ 1 " 1 ^ can be related to the others by 
the AGK cutting rules Q. 



3.1 The four-gluon impact factor 



For the diffractive cross-section, denoted <To above, we will use the result of pO 
where the four-gluon impact factor of the photon is calculated in the asymptotic 
limit. This impact factor is given by multiply cut Feynman diagrams where the 
four gluons couple to a photon via a quark loop (See Fig. f|b). cr is then given by 
convoluting this impact factor with the part of the proton's impact factor which 
contributes to pomeron exchange (That is, with one pomeron-proton vertex on 
each side of the final state cut, the first term of the second row in Fig. ||). The 
proton's impact factor is discussed more detailed in the next section. There, we 
will relate it to the impact factor that is used in single pomeron exchange and to 
the phenomenology of high energy elastic and diffractive scatterings. 

We will now explain why the multiply cut photon impact factor is relevant for the 
calculation of <Jq. The main argument is that the leading part of the amplitude for 
pomeron exchange is purely imaginary. According to the Cutkoski rule, it is then 
given by its s-discontinuity times a factor i/2. For the 0(a 2 s ) amplitude that we 
are considering, taking the s-discontinuity means that we have to take into account 
all planar diagrams (where the gluons do not cross or interact) and put a cut in 
between the two gluons. The practical meaning of the cut is that all particles that 
cross it are on mass shell. When the amplitude is multiplied with its conjugate, we 
will end up with planar diagrams with four-gluons exchanged between the proton 



4 This does, in no way, lead to a full determination of Re„4g^ LLs '. It will give deter- 
mination of the part of Re*4 g ^ LLs \ which contains one gluon from each pomeron, so that 
Re.4 g ^ LLs ^Re„4g[ l NLLs ' ) grows with energy, as fast as the other terms. We will in this way, not 
learn anything of e.g. the vertex corrections for the reggeised gluon. 
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Figure 4: (a) One of the four Feynman diagrams that contribute to the two- 
gluon impact factor of the photon, D( 2; o)(k 1 , k 2 ). (b) One of the 16 contributions 
to the impact factor D( 4 . )(k 1 , k 2 , k 3 , k 4 ). The triple-discontinuity is relevant for 
the diffractive cross-section, since the amplitude for pomeron exchange is purely 
imaginary, to leading order. Note that cut particles are on mass-shell, while the 
gluons (and the photon) can be virtual. 
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and the quark loop and with cutting lines in between all the gluons. The sum of 
the upper parts of these diagrams will then directly correspond to the multiply cut 



photon impact factor defined in EO] 



Before describing the explicit expression for the photonic four-gluon impact factor, 
it is useful to discuss the corresponding one with only two gluons (Fig. |]a). This 
impact factor is given by planar diagrams where two gluons couple to the quark 
loop in all possible ways and with a cut in between the gluons. It has the form of 
a colour tensor, with index a\ and a 2 for gluon one and two respectively, and is a 
function of the transverse momenta (In Breit frame), ki and k 2 , of the gluons. It 
is given by 

D§g(k!, k 2 ) = g 2 8 aia2 l - [ J D (2;0 )(k 1 ) + £> (2; o)(k 2 ) - D (2;0) (ki + k 2 )] (2) 
where D/ 2 -,o) (k) denotes the forward impact factor which, after the loop integration, 



has the analytic form [21 



Dm (k) = ^ £ e f J o da I dy a(1 _ a)Q2 + y{1 _ y)k2 (3) 

Here, Q 2 is the virtuality of the photon, a is the energy fraction of one of the 
quarks and y is a Feynman parameter introduced in order to do the transverse 
momentum integration. The sum is over quark flavours. This expression is valid for 
transverse polarisations of the photon and for massless quarks. For simplicity, we 
will disregard longitudinal polarisation and quark masses. This impact factor has 
the important property that it goes linearly to zero whenever one of its arguments 
goes to zero. The interpretation is that when the gluon becomes soft, it has a long 
wavelength and can not resolve the quark anti-quark pair. Since the quark pair is 
in a singlet system, the gluon does not interact. 

Later on, it will be useful to write the two-gluon impact factor of the proton, in 
the form 

0r 2 (ki,k 2 ) = /^_5 aia2 / 2 ( kl ,k 2 ), (4) 

where N is the number of colours (N = 3). With this definition, we can identify 
the function f 2 in the forward limit, with the gluon distribution in the proton. The 
cross-section for 7*-proton scattering, to 0(a 2 ) is given by 

a {2) = D 2 ®(j) 2 . 

The sub-script specifies that this is a contribution to the one gluon exchange process 
and the (2) denotes the order in a s . The convolution procedure is defined by 
summation over colour indices and integration over the transverse momentum loop, 
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Figure 5: The assumed form for the four-gluon impact factor of the proton. Only 
the first term contributes to the diffractive cross-section, ao. The other two terms 
are needed in the consistency check for the AGK rules (In the Appendix). 



with the measured 

,„ _ d 2 k, d 2 k 2 (27T) 2 <5(k 1 + k 2 ) 

«M = , n „ , n ,o TaU 



(2tt) 2 (2vr) 2 k\k\ 
For the momentum transfer, p±, we can write the differential cross section 

da {2) r 

= irg 4 N J d/C 2 <f(ki - p ± )/ 2 (ki, k 2 )D 2 (k h k 2 ) (5) 

We can now write the explicit form of the four-gluon impact factor [^(J 

I>g.g- a *(k 1> ka J k 8 ,k 4 ) = (6) 

with 

J D^ 0) (ki,k 2 ,k 3 ,k 4 ) = D (2; o)(k 1 + k 2 + k3,k4) + J D( 2 ;o)(k 1 ,k 2 + k3 + k 4 ) 

- D (2 . 0) (ki + k4,k 2 + k3), (7) 
J Dg ;0) (k 1 ,k 2 ,k 3 ,k 4 ) = D (2 .o)(k 1 + k3 + k4,k 2 ) + D (2 . 0) (k 1 + k 2 + k4,k 3 ) 

- D( 2;0 )(ki + k 2 ,k 3 + k 4 ) - D (2;0 )(ki + k 3 ,k 2 + k 4 ). 

Here, g is the strength of the strong coupling and d is a colour tensor defined by 
d abcd = tT ^a t b t c t d^ + tT ^ t d t c t b t a^ where fa denote the SU(3) generators. It will be 

useful for the following, to separate -D( 4; o) into two parts. One part, with the gluons 
Note that the convolution procedure has the dimension of cross-section. 
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one and two (and consequently gluons 3 and 4) in a symmetric colour state, and 
one part for the anti-symmetric colour state. This separation is unique ((+) for 
symmetric and (— ) for anti-symmetric). 

£(4 ia r a4 ( k l> k 2> k 3> k 4) = ^r 2 ^ 4 ^)^ k 2, k 3, ^)+9^ a:ia W[; ) 0) (k 1 , k 2 , k 3 , k 4 ) 

(8) 

with 



^aia 2 a 3 a4 _ ^aia2«3«4 _|_ ^a 2 aia 3 a4^ jc^ 

^01020304 _ (1203014 ^020.1030-4 ^ 

D^^MMM) = ^ ;0 )(ki,k 2 ,k3,k 4 ) + J Dg. 0) (k 1 ,k 2 ,k 3 ,k 4 ) 

^S)( k i, k 2,k3,k 4 ) = D^. 0) (k 1 ,k 2 ,k 3 ,k 4 ) - ^g ;0) (k 1 ,k 2 ,k3,k 4 ) 



The impact factor -D( 4;0 ) vanishes linearly whenever ki or A; 4 becomes soft. This is, 
in general not the case for k 2 and A; 3 . This can be interpreted as that the gluons 1 
and 4 interact with a colour less quark pair, while gluons 2 and 3 interact with the 
quarks in an octet state. One can however check that the symmetric part, -D^q), 
vanishes with any of its argument. Further, it is invariant under the permutation 
of the momenta of the gluons. We know this because it can be written as the sum 
of permutations of _D( 2;0 )(ki, k 2 + k 3 + k 4 ) and — _D( 2;0 )(ki + k 2 ,k 3 + k 4 ). (Note 
that L> (2 . )(ki,k 2 ) = L> (2;0 )(k 2 , ki).) 

The diffractive cross-section a , which will be determined in the following, will 
have all of its contribution given from D^) y In our treatment, this will be the 

case also for the non- diffractive <x 2 and <j[^ , since these are determined from a and 
the AGK rules. This means that we are assuming that the non-diffractive cross- 
sections have additional contributions other than from the multiply cut impact 
factor. This could, for example be from Feynman diagrams with triple gluon 
vertices, where two gluons merge into one before they are absorbed by the quark 
pair. One can speculate that such contributions are important, because they allow 
the "second" exchanged gluon to, in all possible ways, interact with the system of 
the quark, anti-quark and the "first" gluon, which in total is in the colour singlet 
state. 

3.2 The diffractive contribution 

For the part of the proton impact factor, denoted 4 (p), which is used in the 
calculation of <To, we assume the form 

^^(ki, k 2 , k 3 , k 4 ) = g 4 f_ Wa 3 a 4 /4(ki, k 2 ; k 3 , k 4 ) (10) 
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The Kronecker deltas assure colour singlet exchange, and the function / will be 
determined in the next section. For now, it is sufficient to mention that it is 
dimension less, it goes linearly to zero whenever one of its argument goes to zero, 
and it has the following symmetry properties: 

/ 4 (ki,k 2 ;k 3 ,k 4 ) = / 4 (k 2 , k 4 ; k 3 , k 4 ) 
= /4(ki,k 2 ;k 4 ,k 3 ) 
= / 4 (k 3 ,k 4 ;ki,k 2 ). 



The cross-section er , for pomeron exchange, is now given by the convolution of the 
two impact factors, -D( 4; o) and 4>a(v)) 

°0 = ^-D(4;0) ® 4>4,(T) 

where the factor of one fourth is because of the double use of the Cutkoski rules. 

Later on, in section ^ we will be interested in the momentum transfer between 
the proton and the quark pair. Therefore, we write in more detail, the differential 
cross-section for a transverse momentum exchange p± 

da 1 8 iV 

— =- = - • TTO 11 

dp\ 4 y 2 K J 

x J d/C 4 5(ki + k 2 - p ± )/ 4 (ki,k 2 ; k 3 , k 4 )D ( ( ^ ) (ki, k 2 , k 3 , k 4 ) 

d 2 k x d 2 k 2 d 2 k 3 d 2 k 4 (2n) 2 5(k 1 + k 2 + k 3 + k 4 ) 
4 ~ (2tt) 2 (2tt) 2 (2vr) 2 (2vr) 2 k\k 2 2 klkl 

For the colour part of the problem, we have used the result (with d + defined in eq. 

m 

S ab Scddt d = 2C 2 f N 

Where C/ = (A^ 2 — 1)/2N is the colour factor for gluon emission from a quark. 



3.3 Non-diffractive contributions 

We will now make use of the AGK rules to relate Eq. (|TTf) to the processes of non- 
diffractive double gluon exchange and single gluon exchange. It is then important 
to note that kx and k 2 are the internal momenta of the pomeron, up to now, on one 
side of the final-state cut, and k 3 , k 4 are the momenta of the pomeron on the other 
side. The AGK cutting rules relate the contribution do with no cut pomerons to 
the contributions cr\ of one cut and cr 2 of two cut pomerons. We will now identify 
(To with the diffractive cross-section <7o, which we have considered to now, and <7 2 
and <7i, with the cross-section for non-diffractive double gluon exchange <r 2 , and 
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the higher order correction to single gluon exchange a\ , respectively. We assume 
further that the relation er 2 : : a = 2 : —4 : 1 is valid on a more exclusive level 
with fixed values of the transverse momenta fcj. 

For clarity in the following two equations, we will change notations: ki = k a i, 
k 2 = k a2 for the gluons that emerge from the pomeron a, and k 3 = kbi, k 4 = kb 2 , 
for the others. In the case of two cut pomerons, the momentum transfer (The sum 
of the momenta on either side of the cut) can be written as either = k al + k bl 
or p± = k al + k b2 . Because of symmetry properties of the impact factors, the two 
possibilities give the same result, and choosing the first one we get 



da 2 1 8 iV 

= 2'* 9 2 (12) 



x 



J dlC 4 5(Ki + k b i - P±)fa(Ki, k a2 ; k b i, k b 2)-Dj2 } (k a i, k a2 , k b i, k b2 ) 



With dK 4 defined in Eq. ([TT]), / 4 defined in Eq. © and D$ Q) defined in Eqs. © 
and (0). For the case of one cut pomeron, we have any of the possibilities = k; 
but again the choice is irrelevant^ With p^ = k a i we get 



X 



j d/C 4 5(k a i - p_L)/ 4 (k a i, k a2 ; k b i, k^D^^, k a2 , k bl , k b2 ) 



It is important to note that we do not consider all exchanges which are in the 
colour singlet state, as pomeron exchanges. Above, the pomeron was described 
by two gluons in colour singlet. For the contributions with two cut pomerons, we 
have, on each side of the cut, one gluon from each pomeron. These two gluons 
could very well also be in a colour singlet state, but we would not consider the 
process to be pomeron exchange. (This contribution is however colour suppressed 
with a factor 1/(N 2 — 1).) It is not certain that the (partonic) process of the 
exchange of two gluons in a singlet state, as opposed to pomeron exchange, should 
lead to a rapidity gap in the final-state distribution of hadrons, since that would 
correspond to (non-perturbative) colour screening, in the hadronisation process 
(See e.g. P^|). We will leave this question open, but refer to the contribution of 
cr 2 as non-diffractive. 

We will now summarise the main assumptions that was needed to derive the main 
result of this section, that is, Eqs. fll^D and (fT3|). We started by considering the 



3 Note that the number of independent choices for p±, coincides with the AGK weights. 
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process of pomeron exchange. There, we assumed that the amplitude is imaginary 
to leading order and we could therefore use the four-gluon impact factor calculated 
in ||20|| . Further on, we made some assumptions on the form of the proton impact 
factor (The part of it which corresponds to pomeron exchange), including some 
symmetry properties. Finally we assumed that the AGK cutting rules are valid so 
that we could relate the process of pomeron exchange to the processes of double 
gluon exchange and a sub-leading correction to single gluon exchange. The AGK 
rules where assumed to be valid for fixed values of the internal momenta (the 
momenta of the two gluons) in each pomeron, while originally |p.l| , they were 
stated as a relation between the inclusive cross-sections. 

In section [|, we will use the Eqs. (|12|) and (|13|) to make some estimates of observable 
effects of double gluon exchange. In the next section, the impact factor of the 
proton, is discussed and determined on a more quantitative level. 



4 The proton impact factor 



We will start with a discussion of the two-gluon impact factor, 2 , of the proton. 
This is related to the conventional unintegrated structure function or the parton 
distribution of the proton, and there are a large number of calculations and mea- 
surements available for its determination. It would, however, be inconsistent to 
make a direct use of the available parton distributions, since these have, in the 
normal case, been fitted to data, without considering effects from multiple inter- 
action. Instead, we make a simple anzats for 2 , by using the measured parton 
distributions at large momentum fractions, x p ~ 0.1, where corrections from double 
gluon exchange should be small, and by assuming an exponential growth towards 
lower x p - values. 

In the second part of this section, we will determine the four-gluon impact factor, 
04, from the square of the parton distribution, times a normalization factor which is 
determined from two separate measurements in proton-anti-proton scattering. The 
uncertainty in 2 , and consequently in 4 will lead to one of the main uncertainties 
in the final result. In a more careful treatment, it is possible to determine these 
functions more accurately than will be done here. We will instead try to keep 
things as simple as possible, with the aim of just an order of magnitude estimate 
in the final results. 



4.1 The two-gluon impact factor 

With our way of defining the two-gluon impact factor, in Eq. (|^), the function 
/2(ki, k 2 ) is, in the forwards limit, given by the number of gluons, n(k), inside the 
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proton, that an exchanged gluon, with transverse momentum k can resolve. 



/ 2 (k,-k) 



n(k) 



(14) 



The exact meaning of "resolve", in this case, is that we require the transverse 
momentum of the exchanged gluon to be smaller than half of the invariant energy 
of one of the quarks, from the photon side, together with the "resolved" gluon 
in the proton, k 2 < s/4. This leads to a value, x(k 2 ), for the smallest possible 
momentum fraction, x p , that a resolved gluon can have. Apart from k 2 , this limit 
(x) depends on the proton-photon CMS energy (W) and the momentum fraction, 
a of the participating quark from the photon side 



In any case, we will always require that x v > x^j a. We will not try to determine 
a from event to event, but instead, in this case, use the typical value a = 0.5. 

To determine the number of resolved gluons, we will assume a constant contribution 
of n = 2 from the region 0.1 < x p < 1. We will further assume a density of two 
gluons per unit of \og(x p ) at the point x p = 0.1 and from there, that it grows like 
x p A . For the value of the exponent, we will consider the two cases A = 0.3 and 
A = 0.4, respectively. (The gluon distributions, xg\(x), are plotted in Fig. ||) 
This results in the following expression for the number of gluons with x p > x 



where, 6 is the step-function. 

Finally, we have to take into account coherence between the gluons inside the 
proton. When the wavelength of the exchanged gluon is at the order of the proton 
size, it cannot resolve the individual gluons and the effective number of gluons 
will decrease. Therefore, we will assume that the effective number of gluons, goes 
to zero, linearly with k 2 and that the suppression starts around a scale /i 2 . This 
effect is taken into account by multiplying n(x) with a factor k 2 /(/j, 2 + k 2 ). The 
suppression scale is set to /i 2 = 0.1 GeV which is a typical value for the momentum 
transfer in elastic scattering of hadrons (See the discussion in the second part of 
this section). This leads, finally, to 



Xp > 



aW 2 



= x. 




(15) 




(16) 



with n defined in eq. fllBD. 
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Figure 6: The assumed gluon distribution, xg(x), for A = 0.3 (lower solid curve) 
and A = 0.4 (upper solid curve). The dashed curve is the gluon distribution of the 
GWB model. 

This function is plotted in Fig. ^ again for A = 0.3 and A = 0.4 and with 
Q 2 = 10 GeV 2 and xb = 0.0001. We note that the two A-values result in almost 
a factor of two difference in the low-fc 2 region. This will be one of our main 
uncertainties. As mentioned before, it would be inconsistent to use the available 
measured parton distributions, since these have been fitted to data, without any 
account of multiple scattering. Here, we have instead fixed the parton distribution 
at a large momentum fraction x p = 0.1, where saturation effects should be small, 
and then assumed a BFKL like x p - dependence. 

The turnover of n{k 2 ) in Fig. [7|, occurs at the point where x = xb/o>. This happens 
at k 2 = x B W 2 /A = Q 2 /A. 

4.1.1 A comparison with the Golec-Biernat Wiisthoff model 



As stated before, our parton distribution is poorly constrained, because the avail- 
able fits are in general made without any concern to multiple interaction effects. 
There are however several saturation models that have been fitted to the struc- 
ture function, and it is in some cases possible to extract the corresponding par- 
ton distribution from these fits. We will use the model of Golec-Biernat and 



Wiisthoff p3"| , |21 |, which is perhaps the simplest one available and which describes 
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Figure 7: The number of resolved gluons n(k 2 ), for A 
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0.3 and A = 0.4, respec- 



the measured structure function very well, over a large range of Q 2 , in particular 
down to the photo-production region. 

Let us start with a minimal recollection of the GBW model. The description of the 
DIS process is as follows, that the virtual photon splits into a quark-anti-quark pair 
which then propagates through the proton and exchanges one or several gluons. 
The 7*-proton cross-section is, in the GBW model given by 



a n (x,Q 2 )= Jd 2 r\V n (r)\ 2 a(x,r 2 ) 



(17) 



where \1/ is the wave function of the photon (Integrated over the energy fraction 
of the quarks), ir is the polarisation of the photon and r is a vector in impact 
parameter space. The dipole cross-section a is, in this model, set to be 



a(x, r 2 ) = (T I 1 — exp 



4i2g(x), 



with 



Ro(x) 



(19) 



1 /rr\ A /2 
Qo \x J 

The parameters are a , x and A, while Q is set to 1 GeV and is no real parameter. 
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A fit to HERA data resulted in the following values [^|: o"o = 23.03 mb, Xq = 
3.04 • 10" 4 and A = 0.288. 

The validity of the GBW model lies in that it is a simple model which, with a few 
parameters, can describe the structure function over a large range of Q 2 . In fact, 
the authors refer to it as a model which contains saturation, but they only briefly 
mention multiple interactions as a possible explanation for saturation f2"4]j . 

It lies however close at hand to interpret the exponential in the dipole cross-section 
(Eq. [18]) as the poissonian probability of not having any interactions. (The prob- 
ability of not having any randomly scattered partons inside an area r 2 , with the 
average density of partons 1 / (4Rq).) The GBW model is thus a model which incor- 
porates shadowing to all orders in the number of exchanged gluons. An expansion 
of the exponent gives us the corrections that appear when each higher order is 
taken into account. 

A formula, which is useful for extracting the gluon distribution, gives the following 
relation between the dipole cross-section and the (gluonic) unintegrated structure 
function T g {x, k 2 ) (Defined as a distribution in k 2 , as in p4| ) 



cr(x, r) 



4rra s 



rd 2 k 


1 - e irk 


1 k 2 





F g (x,k 2 ). (20) 



This relation can be inverted, but if we use a as it is, then T will not correspond to 
the gluon distribution. That would be true only if we knew that the dipole absorbs 
one and only one gluon in each collision. 

It is however possible to get the gluon distribution of the GBW model by ex- 
panding the dipole cross-section in orders of exchanged gluons, that is, expanding 
the exponent of a . The first term, then would give us the structure function, if 
multiple interactions were absent. This structure function would fit data poorly, 
especially in the low-Q 2 region, but it would instead directly corresponds to the 
gluon distribution. 

Using the first term in an expansion of a we get 

r 2 4VTO, fd 2 k^ -Vi (XjJfe2)> (21) 



CO 



4R 2 (x) 



rd 2 k 


1 - e ir k 


1 k 2 





where the superscript indicates that is calculated from single gluon exchange. 
This equation is inverted to give 

where 6 is Dirac's function. We thus find that the first term in the expansion of a 
gives a distribution peaked at zero transverse momentum. The full a, on the other 
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hand, gives a broader distribution, which is peaked at a finite k 2 . (When xg(x, k 2 ) 
is integrated over all k 2 we get however the same result in both cases.) 

For the gluon distribution, we finally get 

xg(x, Q 2 ) = f F}{x, k 2 )dk 2 = °» I (23) 

The gluon distribution that we find from the GBW model is thus independent of 
Q 2 . This could be expected, since there is no transverse momentum evolution of 
the parton distribution incorporated into the model. The GBW model thus fails 
to describe the structure function in the high Q 2 limit, but the fit shows that the 
strong Q 2 dependence in the lower region can be explained as a saturation effect, 
and is not necessarily related to the usual DGLAP scaling violation mechanism. 

Inserting the fitted parameters we get 

xg(x) ps x . (24) 

a s 

This is plotted in Fig. |B| for a s = 0.3. Note that x in the GBW model corresponds 
to xb of the e-p scattering. The momentum fraction of the interacting gluons is 
not well defined in the model, but is in general somewhat larger than xb- This 
correction would increase the gluon distribution. We conclude that our assumed 
gluon distributions, with a relatively large exponent and independent of Q 2 , are 
supported by data. 



4.2 The four-gluon impact factor 

Now that we have <p2 in the forward limit, we can continue with the discussion of 
how it can be used to determine the two-pomeron-proton form factor, 4 (-p). We 
will rely on the picture that <fi 2 corresponds to a pomeron-proton vertex, and </> 4 (p) 
is, roughly speaking, given by the square of 02- We will also use the fact that in 
high energy elastic and diffractive scattering, one observes low values (well below 
1 GeV 2 ) of the typical momentum transfer. In this article, we will be interested 
in the momentum transfer in DIS processes, and for p 2 ± values above 1 GeV 2 . 
With all the transverse momenta fcj larger than the typical pomeron momentum, 
we can write the function / 4 (ki, k 2 ; k 3 , k 4 ), defined in eq. fllDD, in the following 
approximative way (assuming k* = 0) 

A(pi +r, -pi;p 2 -r, -p 2 ) = / 2 (pi + r, -pi)/ 2 (p2 - r, -p 2 )P(r) 

« n{p 1 )n(p 2 )P{r) (25) 
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Here, we have made the substitutions 



ki = pi + r 

k 2 = - Pl (26) 
k 3 = P2 - r 
k 4 = -P2 

and the function P(r) carries all the dependence on the pomeron momentum and 
a possible normalisation constant. This function will have the effect of an overall 
normalization factor, in the final result 

/fl, T 
p(r W- (2T) 

For the function P(r) we assume the simple exponential form 

P(r) = ae-^ r \ (28) 

which gives 

The value of the parameters a and b will be determined from two different mea- 
surements from high energy proton-anti-proton scattering. 

The first one, we refer to, is a measurement of the exponential slope for the mo- 
mentum transfer, t, in elastic pp scattering. In our model this is described by 
the convolution of the two-pomeron form-factor, (pA(v), with itself. With the ap- 
proximative form in eq. fl25|) , the transverse momentum distribution in pomeron 
exchange is then given by P 2 (r) cx e~ br . Therefore, our parameter b is directly 
given by the exponential slope in pp scattering. 

This slope varies however weakly with the total energy. We will use the re- 
sult |25j from UA4 at CERN, where protons (and anti-protons) were scattered 



at yfs = 546 GeV. This energy is perhaps most relevant for the HERA experi- 
ment. The measured slope actually varies with the value of t, it becomes flatter 
for higher values of —t. Since we want to fix a normalization, the low-t region, 
which dominates, should be more relevant. In [25], a fit is made to the mea- 



sured distribution, and in the region —t < 0.15 GeV 2 , the result for the slope is f\ 
b = 15.2iO.2GeV" 2 . 

The other measurement is that of the so called "effective" cross-section a e s in pp 
scattering. This cross-section can be defined from the cross-section <ji for hav- 
ing two independent high-p^ scatterings (Resulting in four uncorrelated high-p^ 

We are making the approximation r 2 = —t. This is of course a very good estimate. 
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jets.) in a single pp collision. The effective cross-section is then given by a rela- 
tion between cr 2 and the square of the inclusive cross-section for a single high-pj_ 
scattering, according to 

In practice, it turned out to be experimentally much more efficient to measure cr e g 
from events, with one of the four jets, replaced by a photon. This measurement 
was made by CDF f7j, and gave the result cx e ff = 14.5 ± 1.7i 2 !3 mD - 

In our model, the cross-section for one gluon exchange, is given, to leading order, 



iV 2 -l7 (2tt) 2 A; 4 

In this case, we can neglect the higher order contribution, since o\ is the 
inclusive cross-section for one scattering. 

We can now determine the cross-section for double scattering via o-i = 2o"o and 

a ° = \^ ® ^ = T [w=i) x (30) 

/ [/4(pi + r, -pi;p 2 -r, -p 2 )] 2 

where we have made the substitution in eq. fl26|). With the assumption that pi 
and P2 are much larger than r, we can make the approximation in eq. (p5|) 



d 2 pi d 2 p 2 d 2 r 

(2tt) 2 ( Pi + r) 2 p 2 (2vr) 2 (p 2 - r) 2 p 2 (2vr) 2 



g 8 ( N 2 V f 2 2 2 rf 2 Pl d 2 p 2 d 2 r 
a ° = 1 \ N 2 — 1 1 J n ^ Pl ' U ^ ^ 
1 o f „o, , d 2 r 



^y p2 <''W (3l) 



This leads to the simple relation 



Cff 2a 2 A a, JP 2 (r)^ 



and it gives us the second relation, for the determination of a 



£ = = £ (33) 



8 We have simplified the notation considerably by suppressing dependencies on the longitudinal 
variables. The final result is not affected by this. 
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From Eqs. (^) and (P^), and with the values b ~ 15 GeV 2 and a c fj « 15 mb ~ 
37 GeV -2 we get a ~ 2.3, and the normalization factor 

N = -. ] « 0.024GeV 2 
VvrOcr eff 

Now, we have everything we need for the calculation of the relative effect of dou- 
ble gluon exchange, on the distribution of the momentum transfer in 7*-proton 
scattering. 



5 Numerical results for the momentum-transfer 

In this section, we will present some numerical results from the model for double 
gluon exchange, that we have developed. We will consider 7*-proton scattering, 
and make an estimate of the relative effects from double gluon exchange, on the 
distribution of transverse momentum, that is exchanged between the proton and 
the quark pair that originates from the virtual photon. 

More specifically, we will make an estimate of the function 



da2 da[ 4 ^ \ I da^ 
dp± dp± J I dp\ 



A(i4)= l^ + rzH l-dr (34) 



We will now write the formulas for the three differential cross-sections in detail. 
For the 0(a 2 ) contribution to one gluon exchange, we have, from Eqs. (|5|) and 

© ^ (2) 

= 4Nna 2 s D 2 (p ± )±n(p ± ) (35) 
ap ± p ± 

For the contribution from double gluon exchange, we get from eq. (0) and the 
substitution in eq. (|26|) 

d(To ir s 1 r d 2 pi d" 2 p 2 r9 „. s , 

9 N 7^ / 7 r^W ^r^d 2 rS( Pl + p 2 - p ± ) x (36) 



dpi ~ 4 y (2nfJ ( Pl + r) 2 p 2 (p 2 - r) 2 p 2 

>(+) 
(4;0) 

,M rd 2 Pl P \_ 



x h(Pi + r > -Pil P2 - r, -p 2 ).D, ( +jL(pi + r, -pi, p 2 - r, -p 2 ) 



« ANira*— / — r - -n (pi n (p ± - p : x 

Pi Pi (Pi - Pi) 

x ^o)(Pi>-Pi>P-L-Pi>Pi _ P-l) 



where we have made the approximation in eq. (|25|) and used the definition in eq. 
(E7|). The normalization constant has the value J\f ~ 0.024 GeV 2 . In the same 
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way, we can write for the 0{af) correction to the single gluon process^ eq. (P~5| ) 



d(?F re 8 1 r cPpi d 2 p 2 2 

-n9 N 7^Z / 7Z I .N9„2 7Z 3?l2 rf r5 (Pi - P-l) x ( 37 ) 



dpi " 2 y (2vr)6 7 ( Pl + r)2p2( p2 _ r )2 p 2 

)(+) 
(4;0) 

i A/" /-rf 2 p 2 



x h(Pi + r, -pi; p 2 - r, -p 2 )D f ( +L(p 1 + r, - Pl , p 2 - r, -p 2 ) 



-8A%a 4 — f —^ L n(p ± )n(p 2 ) D ( ^l(p ± , -p ± ,p 2 , -p 2 ) 

P± J P2 



With the relation 



D$o) (Pi , -Pi , P2, -P2) = 2 J D (2; o) (pi) +2D( 2; o) (P2) ~D(2fl) (pi +P2) -D {2 . 0) (pi -p 2 ) 

derived from the Eqs. @, (0) and @, we can, finally, write the result for the 
correction function, defined in eq. ([34]) 



2 A/" 



<i 2 pi n(pi] 



P± n(p± - pi, 



x 



X 

-2 



Pl £(2;0)(p±) [n(p±) (p± -Pl) 4 
2-D(2;0)(Pl) + 2D (2;0 )(p_L - Pl) - £>(2;0)(P±) ~ ^(2;0)(2pi - P_L 

2-D( 2; o)(Pi) + 2D( 2 ;o)(p±) - -D( 2; o)(p± + Pi) - -D(2;0)(P± - Pi)]} • 



(3? 



In Fig. ||, two different estimates of the correction function A(pj_) are plotted for 
electron-proton DIS at xb = 0.0001 and Q 2 = 10 GeV 2 . The integration over pi 
was done numerically. For the two-gluon impact factor -D( 2; o) (See eq. (§)) the y 
integration was done analytically, while, for the a dependence, the constant value 
a(l — a) = 1/6 was inserted^. The larger estimate of A in Fig. || is given by 
A = 0.4 and a constant coupling a s = 0.3. The smaller estimate is from the choices 
A = 0.3 and a s = 0.2. 

The lower choice for the coupling, corresponds to a scale of the order 20 GeV 2 . 
Though we are interested in the distribution at rather large p± values, it seems 
likely that the correction to the p^-distribution is mainly due to an additional 
rather soft gluon, with a transverse momentum of the order 1 GeV 2 . In case, this 
is the correct scale as the argument of a s for the second gluon, the larger value for 
the coupling, would be more appropriate. 

We see in fig. (|5p that the correction function has its maximum at around 15 GeV 2 
with the maximal corrections of 6 % and 23 %, respectively, in the two estimates. 
Towards higher p± values, the correction decreases, but continues to be positive. 

9 This is by no means all of the 0(a 4 ) corrections to the one gluon exchange process (Not even, 
in the asymptotic limit). The idea is that, it is all of the 0{a s ) correction, in the asymptotic 
limit, which is proportional to n 2 . 

10 The average value of a(l — a) in the range [0 : 1]. This is, of course, an approximation. 
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Figure 8: Upper and lower estimates for the correction function A(p 2 ± ), (Eq. (|3~4D). 



For the low p± values, there is a large negative correction. In this region, the 
correction from three-gluon exchange should be positive. In the case that the 
negative corrections in this region are large (At the order of one, or greater) they 
could not be trusted since the positive correction from three gluon exchange would 
also be large. When it comes to the negative corrections in the low-p^ region, one 
is, for any treatment with a fixed number of exchanged gluons, always dependent 
on some kind of lower cut-off. 

These negative corrections, at the order of one, are however, a signal for an in- 
teresting phenomenon that have been expected from other models (See e.g. If2q|), 
namely, saturation in the transverse momentum distribution. If multiple gluon ex- 
change to all orders, are taken into account, one expects that the, p_i_-distribution, 
which is peaked at zero for one gluon exchange, will instead become flat below 
a lower limit, A. One can make the interpretation, that the density of gluons is 
so large, in this region, so that in almost every scattering event, at least one of 
them is absorbed by the quarks. The limit, A, is expected to grow, towards lower 
a^-values [^6 . 



We will now list the main approximations that have been used and that can affect 
our result for the p± distribution. The severity of each of these approximations 
is estimated by the factor, p, that the (positive) maximal value of the correction 
function, A is changed, if we would do the approximation, in another, but still 
reasonable way. In this respect, we will refer the lower estimate as the "Default" 



26 



choice (See Table [I].). 

There are three entities that are directly proportional to the final result. These 
are the normalisation factor A/", the square of the coupling constant a 2 and the 
number of resolved partons n(k). 

For the normalisation factor, A/", the experimental uncertainty in cr e ff of about 10 %, 
will result in the same relative error in A. An error, perhaps equally important, 
should arise from our simple physical assumptions in the determination of M . 
Further, the two choices of a s = 0.2 and a s = 0.3 give a relative factor p = 2.25. 

The discussion of errors on the number of partons n(k) is more involved, but 
important. From Eqs. (|i~5|) and (|16"D, one can see that there are three parameters 
that determine n(k). These are the exponent A, the gluon density at x v = 0.1 
(This was set to 2) and the scale at which, coherence effects in the proton sets in, 
p 2 = 0.1 GeV 2 . 

The scale p? has rather small effects on the maximal (positive) value of A. An 
increase to the value p 2 = 0.5 GeV 2 , reduced the maximal positive correction from 
~ 6% to ~ 4% (p = 0.67). The maximum was shifted to larger p± values, with 
approximately 5 GeV 2 . 

For the dependence on A, by changing the value from 0.3 to 0.4, the number of 
resolved gluons, in the low- A; 2 region is more than doubled (See Fig. [/]). This 
results in an doubling of the maximal correction, that is p = 2. 

There is another uncertainty on our gluon density n(k), which is related to the 
fact that we have used a formalism which is valid in the asymptotic s — > oo limit. 
In the derivation of n(k), we introduced the lower limit x(k 2 ), for the momentum 
fraction, that a resolved gluon could have. This limit was derived from a local 
invariant energy squared s and we required, for the exchanged gluon that k 2 < s/4. 
One should expect that, compared to the asymptotic result, there is a suppression 
also, in a range below this value (In the case of time-like cascades, see e.g. f27j). To 
estimate the significance of this approximation, the function n(k) was reevaluated 
with the limit k 2 < s/8. This lead to a reduction of n(k) of about 20 %, over the 
whole range, except for p± values below 2 GeV 2 , where n was not changed. In this 
region, the limit x is insignificant because it becomes smaller than xb- The effect 
on the correction function A, was that it had the same shape and maximal value, 
but was shifted to lower p± values, a distance of 5 GeV 2 . So, in this case, we had 
p=l. 

It might be interesting to check how the maximal value of the correction, A, varies 
with the values of xb and Q 2 . Therefore, separate runs were made for Q 2 ) = 
(10- 3 , 10 GeV 2 ), (10- 5 , 10 GeV 2 ), (lO" 4 , 20 GeV 2 ) and (10~ 4 , 5 GeV 2 ) (See Table 
HI for the results for p). The conclusion is that the effect of double gluon exchange 
increases rapidly towards, lower values of both xb and Q 2 . Another observation 
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from these runs is that the maximum always occured at p\ values slightly above 
Q 2 . The Q 2 -dependence of the maximal correction, seems to be consistent with the 
expectation that multiple interactions is a higher twist effect. A l/<5 2 -behaviour 
could also be explained by the fact that A is dimension less, but proportional to 
M which has the dimension [energy] 2 . (The dependence on /x 2 was observed above 
to be weak, perhaps logarithmic). 



Entity 


Alternative procedure (Default) 


P 


M 


[Experimental errors] 


±10 % 


M 


[Theory assumptions] 


±10 % 


a s 


Set to 0.3 (0.2) 


2.25 


n 


^ = 0.5 GeV 2 (0.1 GeV 2 ) 


0.67 


n 


A = 0.4 (A = 0.3) 


2 


n 


x determined by k 2 <s/8 (k 2 <s/4) 


1 


x B 


io~ 5 (io- 4 ) 


1.9 


x B 


io- 3 (10~ 4 ) 


0.5 


Q 2 


5 GeV 2 (10 GeV 2 ) 


1.8 


Q 2 


20 GeV 2 (10 GeV 2 ) 


0.5 



Table 1: Estimate of various uncertainties in the final answer, p is the relative 
effect on the maximal correction. In the second part of the table, the dependence 
on and Q 2 , is displayed. 

There have been several measurements and analyses of jet-distributions at HERA, 
in the region of low xb and Q 2 > 1 GeV , that is relevant for the discussion here 
(E.g. [p8|l ). The high-p_i_ jet activity, has been described by the definition of jet- 
rates: The fraction of the cross-section, with a number of jets above some high 
p_i_-value. There has also been specific studies of the fraction of the photon energy, 
that is carried in these jets. In this way, one determines if it is more likely that the 
jet originates from a quark (high energy fraction) or a gluon, which is more likely 
to be in the proton direction (low fraction of the photon energy). 

The result, for the momentum transfer, that has been presented here, can not be 
directly compared with the measured jet-distributions. The main reason is that the 
transverse momenta of the quarks was integrated out. It is, of course, not necessary 
to do this, but the ^-distribution of the individual quarks, in the high-p^ region, 
is very much dependent on the energy that is available for the quark-anti-quark 
system. In the leading log(s) formalism, that we have used, there is no information 
on this energy (It is assumed to be infinite), and one would have to determine it, 
by other argument sQ. 

n This means that we have overestimated the cross-sections, in the treatment above, but this 
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So far, there has not been any direct measurements of the momentum transfer to 
the quark pair. In the case of very high momentum transfers, p\ > Q 2 , one could, 
however, be rather certain, that this will lead to two jets, with transverse momenta 
~ p±, originating from one of the quarks, and the gluon which takes the recoil (The 
largest recoil, in the case of double gluon exchange). The other quark, with a lower 
transverse momentum, should then be observable as the photon remnant. Our 
correction to the momentum transfer, would therefore lead to similar corrections 
(enhancements) in the so called "resolved-photon" events, with high p±. 

Another interesting observation, that one could perhaps make at HERA, is re- 
lated to the negative correction, in the region of low momentum transfers. If the 
saturation limit (A above) is sufficiently large, it should be observable in the p±- 
distribution of the quark jets, in a frame were the proton and the photon have zero 
transverse momentum. This would be possible, only if A is large enough, in the re- 
gion of the lowest x^-values, available at HERA, so that the effect is not drowned 
by the hadronisation process. From the larger estimate, above, this saturation 
should be expected to occur at A 2 ~ 2 GeV 2 . 



6 Summary 

We have studied the process of double gluon exchange, in non-diffractive 7*- 
proton DIS. A new model was developed, based on previous descriptions of diffrac- 
tive events, as the perturbative exchange of two gluons, originating from one 
pomeron [HJ. The relation to the non-diffractive process, was determined from 



a general result from Regge theory: The AGK cutting rules [IT]. The validity of 



this procedure has, however, not been fully proven (In the leading log(s) expansion, 
these contributions are next-to-next-to-leading.). 

As a first application of the model, we studied the corrections to the transverse 
momentum transfer to the quark- ant i- quark pair, at the photon vertex. We found 
significant enhancements (At the order of 10%) of the cross-section for large mo- 
mentum transfers, and large negative corrections, (~ 1) for momentum transfers 
close to zero. The prediction for the negative corrections is however not reliable, 
but signals saturation in the transverse momentum distribution, in an all orders 
(in the number of exchanged gluons), treatment. The possibility of the observation 
of these effects in the measured jet-distributions at HERA, was discussed. 

The results presented here, should be taken as order of magnitude estimates, and a 
description of the general behaviour. It should be possible to make the model more 
accurate, in order to make direct comparisons with data. Some of the necessary 



error is small inclusively over the quark momenta and should rather well factor out in the relative 
correction A. 
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improvements, like a new fit of the parton distributions, are straight forward, but 
involved. The most effective way of doing this would probably be to include the 
perturbative multiple interaction as a new feature in existing models for DIS. 
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A consistency check 

In section [J[ we could make use of the Cutkoski rule since the amplitude for 
pomeron exchange is purely imaginary to leading order. This is expected for am- 
plitudes which describe the exchange of a symmetric colour state. We will now, 
with the same method as o"o was calculated, calculate the contribution to ct 2 which 
is due to just the imaginary part of the amplitude for two gluon exchange. For 
the exchange of symmetric colour states, this is the leading contribution. Conse- 
quently, our method should give the correct AGK relation <r 2 : <Jq — 2 : 1, for the 
parts of o~ 2 which are due to the exchange the symmetric colour states and the cor- 
responding part of do- The main argument is that both <jq and the symmetric part 
of o"2, have their contributions given from the symmetric part, D^(ki, k 2 , k3, k 4 ) 
(Defined in Eqs. (||) and (0)), of the photon impact factor, and that this function 
is invariant under permutations of its arguments. 

To do this, we need the full impact factor of the proton (Fig. [5]), whereas, before, 
we used only the part with one pomeron-proton vertex on each side of the final- 
state cut. We choose a form which is explicitely symmetric in the simultaneous 
exchange of the colour and momentum of any two gluons. 

0r 2a3a4 (ki,k 2 ,k 3 ,k 4 ) = (39) 
0™ a4 (k 1? k 2 , k 3 , k 4 ) + 0^r a4 ( kl ' k 3, k 2 , k 4 ) + ^ a3a2 (k 1; k 4 , k 3 , k 2 ) 

with 4 (-p) defined in eq. (|T0D. 

When the impact factor above is convoluted with the photon's impact factor, we 
will get three terms, arising from the three terms in 4 . The first of these is what 
was discussed before, the cross-section cxo, while the sum of the other terms is the 
contribution to cr 2 . As discussed before, for the contributions to cr 2 which are due 
to the exchange of symmetric colour states, we can again make use of the multiply 
cut photon impact factor, but now, convolute it with the second and third terms 
of 4>4. To do this, we have to make projections into the symmetric colour states. 
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In the case of cx 2 , we can do this by putting in a symmetry projector on the colours 
of gluon one and two. This part of the impact factor has already been singled out 
in eq. (§). We can now see that, in the convolution with the second or third terms 
of 04, we get the same colour factor: 



/ 



d a + bcd 5 ac 5 bd = df d S ad S bc = (c f - ^) NC 

Further, the two terms have the factors 

J D(+)(k 1 ,k 2 ,k3,k 4 )/(k 1 ,k 3 ;k 2 ,k 4 ) 

and 

J D(+)(k 1 ,k 2 ,k3,k4)/(k 1 ,k 4 ;k3,k 2 ) 

respectively. Since is invariant under permutations of the momenta, these 
two factors are equal to each other (After renumbering) and also equal to the 
corresponding factor which appears for <jq. 

The corresponding colour factor for gq is 

dt%bS cd = INC) 

From the AGK rule, this would mean that the total contribution to <r 2 should have 
the factor 4iVC|. Above, we have calculated the contribution from the exchange of 
symmetric colour states. Consequently, the validity of the AGK rule, gives for the 
anti-symmetric exchange part of <r 2 , that it should come with the same momentum 
dependence, and with the colour factor 

ANC 2 f -2(c f - —] NC f = N 2 C f 



f — Z O f iV O f 

} \ t 2NJ i 

We will now calculate the colour factor for the part of <jq which corresponds to 
the symmetric colour exchange in <r 2 0- For o"o, the colour symmetry projector 
should be connected with one gluon from each pomeron, e.g. the gluons 2 and 3 
or 1 and 3. There is also a second requirement, namely that the two colour indices 
that are symmetrised over are directly connected to each other through one quark 
line. This is the case for the symmetrisation procedure that was used for the cr 2 
contribution: 

( jaia 2 a 3 a4 ^ ( jaia 2 a :i a4 _ \ ^ ( jaia 2 a 3 a 4 _|_ ^a2«i«3«4^ 

The symmetrisation is done over neighbouring cites of the tensor d. 
For o"o, this can be done in the following way for the first term in Eq. 

^a\a 2 aza,4 ^ \ ^aia2U3i4 _|_ ^a,\aza 2 a^ 



12 It is acknowledged that the exact definition of this procedure is perhaps not known. 
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and for the second term 

( ja 2 a 1 a 3 a4 ^ \_ ^a 2 aia 3 a 4 _|_ ^02030104^ 

2 

After this procedure, and multiplying with S aia2 S a3a4 , we get the same colour factor, 
(C f - l/2N)NC f , also for the a contribution. 

Thus, after the symmetrisations, we get the same contribution to a from the first 
term of 04, as each of the contributions from the other two terms to a^. That 
is, the symmetrised contributions satisfy the AGK relation a : a 2 = 1 : 2. This 
result depends on the form of 04 that has been used, but we have shown that it is 
possible to choose a form, inspired by an underlying picture of two pomerons that 
couple in all possible ways, so that the AGK rule is satisfied. 
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